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We present a study of excited charm-strange baryon states produced in e + e~ annihilations at 
or near a center-of-mass energy of 10.58 GeV, in a data sample with an integrated luminosity of 
384 fb -1 recorded with the BABAR detector at the PEP-II e + e~ storage rings at the Stanford Linear 
Accelerator Center. We study strong decays of charm-strange baryons to A+ K%, A^K" , A+A"~7r + , 
h.t K®n~ , ' KgTv~ tt + , and A J K~ n + tv~ . This study confirms the existence of the states E c (2980) + , 
H c (3077)+, and H c (3077)°, with a more accurate determination of the H c (2980) mass and width. 
We also present evidence for two new states, H c (3055) + and H c (3123) + , decaying through the 
intermediate resonant modes E c (2455) ++ if~ and E c (2520) ++ if ~ , respectively. For each of these 
baryons, we measure the yield in each final state, determine the statistical significance, and calculate 
the product of the production cross-section and branching fractions. We also measure the masses 
and widths of these excited charm-strange baryons. 
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I. INTRODUCTION 

With the observation of the f2* (ess) baryon [1] by 
the BABAR Collaboration, every predicted SU(4) ground- 
state single-charm baryon has been experimentally ob- 
served. Several excited A c (cqq), S c (cqq), and S c (esq) 
baryons have also been experimentally observed [2] . The 
spins and parities of these are assigned based on a com- 
parison of the measured masses and natural widths with 
predictions of theoretical models. 

Both the BABAR and Belle collaborations have searched 
for ground-state double-charm baryons decaying to the 
final state A+A'~7r + [3 5]. These searches reveal no ev- 
idence for such states. However, the Belle Collabora- 
tion finds evidence for two excited charm-strange baryon 
states, S c (2980) + '° and S c (3077)+'°, decaying strongly to 
A^K®ir~ and K^K~tt + [5]. Although these new states 
have the same or similar decay modes as those used in the 
search for weak decays of double-charm baryons, they are 
identified as charm-strange states based on the measured 
masses, natural widths (which indicate strong decays), 
and charges of the members of the isospin doublet. 

Previously known excited S c baryons have been ob- 
served only in decays to a lower-mass S c baryon 
plus a pion or photon. In contrast, the S c (2980) +, ° 
and S c (3077) +, ° baryons are observed in decays in 
which the charm and strange quarks are in separate 
hadrons. The observed decay modes may have impli- 
cations for the internal quark dynamics of these new 
states. Several excited charm-strange baryons with J p = 
{(l/2) ± , (3/2) ± } are predicted, with masses ranging from 
about 2800MeV/c 2 to 3150MeV/c 2 [6, 7]. Other au- 
thors [8-10] consider J p = (5/2) + states and radial exci- 
tations, and use the measured natural widths and decay 
modes in assigning possible quantum numbers for the 
new 5 c (2980)+-° and S c (3077)+'° states. 

In this paper, we report measurements of excited 
charm-strange baryon properties based on a data sam- 
ple corresponding to an integrated luminosity of 384 fb -1 
recorded with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage rings at the Stanford 
Linear Accelerator Center. We search for decays to 
the three-body final states A+AT - 7r + and A+AT!?7r _ , the 
two-body final states A+if° and A+A' - , and the four- 
body final states h+K®n~ir + and A+ K~ir~ir + . Signifi- 
cant signals are found only in decays to three-body final 
states. We confirm the existence of the states S c (2980) + , 
S c (3077) + , and S c (3077)°, with an improvement over ex- 
isting measurements of the S c (2980) + mass and width. 
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For the three-body final states, we also search for 
decays through intermediate resonant £ C (2455) ++ A'~, 
£ C (2455)°A1!, E c (2520) ++ AT-, and £ C (2520)°A:° chan- 
nels. The £ c (2455) and £ c (2520) baryons decay ex- 
clusively to A+7r. We find evidence for two additional 
new states S c (3055) + and S c (3123) + decaying through 
the intermediate-resonant channels £ C (2455) ++ AT _ and 
£ C (2520) ++ AT _ , respectively. We measure the mass, nat- 
ural width, yield, and the product of production cross 
section and decay branching fractions when there is evi- 
dence for an excited charm-strange baryon. Also, where 
applicable, the intermediate resonant decay fractions are 
measured. 



II. ANALYSIS OVERVIEW 

This analysis relies primarily on the charged-particlc 
tracking and particle-identification capabilities of the 
BABAR detector. A detailed description of the BABAR 
detector is presented in Ref. [11]. The charged-particle 
tracking system consists of a five-layer double-sided sili- 
con vertex tracker (SVT) and a 40-layer drift chamber 
(DCH). Discrimination between charged pions, kaons, 
and protons relies on ionization energy loss (dE/dx) in 
the DCH and SVT, and on Chcrenkov photons detected 
in a ring-imaging detector (DIRC). A CsI(Tl) crystal 
calorimeter is used to identify electrons and photons. 
These four detector subsystems are mounted inside a 1.5- 
T solenoidal superconducting magnet. The instrumented 
flux return for the solenoidal magnet provides muon iden- 
tification. 

We produce samples of simulated events using the 
Monte Carlo (MC) generators JETSET74 [12] and EVT- 
GEN [13] with a full detector simulation based on 
GEANT4 [14] . We produce about four million simulated 
e + e~ — > cc events in which at least one of the primary 
charm quarks hadronizes into an excited charm-strange 
baryon that decays according to one of the studied decay 
channels. All particle decays are generated according to 
phase space. Reconstruction efficiencies are estimated 
based on the excited charm-strange baryon kinematic 
distributions from JETSET74. The simulated samples 
are also used to estimate measured invariant-mass reso- 
lutions. 

In data, we study the invariant-mass region between 
2.91GeV/c 2 and 3.15GeV/c 2 for potential S c states. We 
do not examine candidates in this region while optimiz- 
ing the selection criteria in order to minimize poten- 
tial experimenters' bias. Selection criteria are chosen to 
maximize the expected significance of signal events aver- 
aged over five reconstructed A+ decay modes: pK~Tr + , 
pAT°, pK®ir~ir + , A-7T+, and Air + ir~ir + . The expected 
significance of each decay mode is estimated based on 
the efficiency measured with samples of simulated signal 
events, and the number of candidates in upper (above 
3.15 GcV/c 2 ) and lower (below 2.91 GcV/c 2 ) E c invariant- 
mass side-bands in data. 
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The selection criteria are based on proton, kaon, and 
pion particle identification, among several other recon- 
struction parameters. Reconstructed tracks in the en- 
tire candidate decay chain arc simultaneously fit with 
vertex constraints and with the A+, A, and K® candi- 
date masses constrained to their world-average values [2] . 
The total x 2 probability of the fitted track vertices in 
the decay chain is required to be greater than 1%. The 
vertex-constrained A+ mass (with no mass constraint) 
is required to be within 1.75 times the mass resolution 
(~ 5MeV/c 2 ) of the world-average A+ mass [2]. The 
vertex-constrained A and K® masses are required to be 
within 6.5MeV/c 2 and fO.OMeV/c 2 , respectively, of their 
world-average values [2] . The measured A and flight 
lengths divided by their errors are required to be greater 
than 2.5 and 3.0, respectively. The momentum of the 
S c candidate in the e + e _ center-of-mass frame (p*) is 
required to be above 2.9GeV/c. This p* requirement sig- 
nificantly reduces background from random combinations 
of tracks. 

For each S c decay mode, the selected candidates in 
data are divided into five samples based on the recon- 
structed A+ decay mode. The signal yields in each sam- 
ple are extracted with an extended, unbinned, maximum- 
likelihood fit. All five samples arc fit simultaneously. The 
signal probability density functions (PDF) share the val- 
ues of parameters that describe the signal shape and, for 
three-body final states, the intermediate-resonant decay 
fractions. For A+ decay mode m, the corresponding por- 
tion of the likelihood function L m has the form 



L, 



n 



(1) 



where is the number of candidates in signal- 

component T), -B^j is the number of candidates in 
background-component k, a and b are the shape param- 
eters for the N s signal-PDF components (Pj 7 ) and the 
iVj, background-PDF components (P b K ), and x l m repre- 
sents the measured masses for the N m candidates in de- 
cay mode m. The mass, natural width, and resonant 
fraction results listed in Table I are the measured shape 
parameters a for the respective signal PDFs. Combining 
the five L m , the full extended likelihood function C has 
the form 



C — cxp 



5 

E 



n 



(2) 

Each parameter is allowed to vary over a range that is 
large enough so that it does not constrain any of the final 
results; in particular, the number of signal candidates is 
allowed to be negative. For each signal-PDF component 
(PJ 7 ), the numbers of measured candidates with A+ de- 
cay mode m (S 1 ^) are summed together to determine the 
total yields listed in Table I or to determine the prod- 
ucts of cross-sections and branching fractions listed in 
Tables III and IV. 



Statistically significant signals for excited charm- 
strange baryons are observed only in the analysis of three- 
body final states. For the two-body and four-body final 
states, we only search for the baryons observed in the 
three-body final states. The S c masses and widths mea- 
sured with the three-body candidates are used as Gaus- 
sian constraints in the fits to samples of two-body and 
four-body candidates by multiplying the likelihood func- 
tion in Eq. (2) by 



exp 



(M /i) 2 (r - 7 ) 5 



2a 2 



(3) 



where M. is the mass and T is the natural width of the S c 
state measured from the fit to three-body candidates, um 
and <7r are the corresponding uncertainties, and fi and 7 
are the constrained signal mean and width parameters. 



III. DECAYS TO A+K--K+ AND A+J^Tr - 

The studies of the three-body final states A+K~n + 
and A+K®ir~ are based on fits to the two-dimensional 
invariant-mass distributions A/(A+/\~7r + ) versus 
A/(A+7T+), and M{A+K s -k~) versus M(A+7r _ ). The 
experimental resolution for these two- and three-body 
invariant masses varies from about 1.0MeV/c 2 to 
2.5MeV/c 2 . With these two-dimensional mass distribu- 
tions, we can incorporate the intermediate resonances 
£ c (2455)++'° and £ c (2520)++'° in the fits. No other 
known intermediate resonances are kinematically al- 
lowed. We fit a region defined by the kinematically 
allowed mass thresholds, and A+K®ir~ and A^K~tt + 
invariant masses up to 3150 MeV/c 2 . Two-dimensional 
mass plots for candidates in the fit region are shown in 
Fig. 1(a) for A+K~tt + candidates and in Fig. 2(a) for 
A+AT!?7r~ candidates. 

For the two-dimensional maximum-likelihood fit we 
use a PDF with background components (P&) and signal 
components (P s ) used to describe four candidate cate- 
gories: resonant combinatoric background, non-resonant 
combinatoric background, resonant signal, and non- 
resonant signal. We use a double- Voigtian resonance 
shape (a non-rclativistic Breit-Wigner distribution con- 
volved with a resolution function that is the sum of two 
Gaussian functions) to describe each peak in invariant 
mass. These peaks include the 5 C states as well as the S c 
states associated with both 5 C decays and backgrounds. 
The double- Voigtian resonance shape does not account 
for possible interference effects due to overlapping res- 
onances. All resolution parameters are fixed to values 
measured in the MC samples. The masses and widths 
of the £ c (2455) ++ and E c (2520) ++ are free parameters 
in the fits, while the masses and widths of the S c (2455)° 
and S c (2520)° are fixed to the world average values [2] 
because of the limited statistical sensitivity of the data 
to these states. The double- Voigtian resonance shape 
that accounts for excited charm-strange baryon decays 
through an intermediate S c resonance is multiplied by 
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a two-body phase-space factor, and the double- Voigtian 
resonance shape that accounts for direct three-body de- 
cays is multiplied by a three-body phase-space factor. 
The non-E c backgrounds are modeled with a threshold 
function proportional to 



M 



— -1 



exp 



«l f I 



(4) 



where M is the mass variable in the distribution of which 
there is a minimum kinematic threshold, t is the value of 
the threshold, and a and (3 are shape parameters. The 
free parameters (a) of the signal-PDF components in- 
clude the masses, natural widths, and resonant fractions 
of the excited charm-strange baryon signals. 

A. A+K--K+ Results 

Projections of the A~^ K~n + invariant-mass distribu- 
tion, for all five reconstructed A+ decay modes com- 
bined, are shown in Figs. 1(b) and (c) with the combined 
components of the fit shown. The fit allows for four ex- 
cited charm-strange baryon states; the masses and widths 
are free parameters in the fit. The data and fit projec- 
tions correspond to M(A+7r + ) ranges within 3.0 natu- 
ral widths of the £ c (2455)++ mass (Fig. 1(b)) and 2.0 
natural widths of the £ c (2520) ++ mass (Fig. 1(c)) [2]. 
These M(A+7r + ) ranges are delineated by the horizon- 
tal lines in Fig. 1(a). A two-dimensional (M(A+7r + ) 

o is 



and M(A+A'~7r + )) binned % 2 probability of 
calculated for the fit using roughly 300 equally sized 
bins each containing over 10 candidates. Four ex- 
cited charm-strange states are observed at the approx- 
imate M(A+K~tt + ) invariant-masses of 2970MeV/c 2 , 
3055MeV/c 2 , 3077 MeV/c 2 , and 3123MeV/c 2 . We use the 
change in maximum ln-likclihood value when a signal is 
removed from the fit, and the number of free parameters 
in each signal PDF component, to estimate statistical sig- 
nificances of > 9. Oct (Aln£ = 81), 6Aa (AhxC = 31), 
> 9.0cr (A In/: = 91), and 3.6 a (Aln£ = 15) for each 
state, in order of increasing mass. These estimated sig- 
nificances do not statistically account for fluctuations 
in background level across the entire range of invariant 
masses over which the study was conducted. This means 
that these significances do not correspond to those of a 
search. Rather, they correspond to significances of sig- 
nals incorporated into the PDF after evidence for the 
corresponding states had been observed. 

The S c (2980)+ PDF includes a nonrcsonant A+K~tt + 
component and a £ C (2455) ++ A'~ component. (The 
£ C (2520) ++ AT~ final state is not kinematically allowed 
for the S c (2980)+.) For the S c (2980)+, (55 ± 7 ± 13)% 
of the signal is found to be due to decays through 
£ C (2455) ++ A'~. (The uncertainties are statistical and 
systematic. The estimation of systematic errors is dis- 
cussed below.) The S c (3077)+ PDF includes a non- 
resonant A+AT~7r + component, a £ C (2455) ++ AT _ com- 
ponent, and a £ C (2520) ++ AT _ component. For the 



S c (3077)+, (95 ± 14 ± 6)% of the signal is found to be 
due to decays through intermediate resonances. A lower 
limit of 80% for this fraction is determined, at 90% confi- 
dence level, from a numerical integration of the posterior 
probability density calculated from the likelihood distri- 
bution and a zero prior for fractions above 100%. Of 
these intermediate-resonant decays, (45 ± 5 ± 5)% are 
through £ C (2455) ++ A'~. Separate fits to the data that 
include A+AT-7T+, S C (2455) ++ AT-, and £ C (2520)++A"- 
PDF components for S c (3055) + and S c (3123) + signal 
indicate that these states only have signals for decays 
through the intermediate-resonant states £ C (2455) ++ A' _ 
and £ C (2520) ++ AT~, respectively. For all measurements, 
the S c (3055)+ and S c (3123)+ PDFs include only these 
intermediate-resonant components. 

The newly identified S c (3055)+ and S c (3123)+ 
baryons are found to have statistically significant sig- 
nals with widths that are larger in value than our mass 
resolution, p* values distributed toward higher momenta 
than the combinatoric background, and consistent signal 
decay rates between the five A+ final states. Four addi- 
tional samples are studied for possible sources of peaking 
background that could be mistaken for any of the ex- 
cited charm-strange baryon signals. These four samples 
are a "wrong-sign" A+AT + 7r _ data sample, a A+ mass- 
sideband data sample, a MC sample of e + e~ — > cc events, 
and a data sample in which the A+AT _ 7r + invariant mass 
is recalculated substituting a pion mass for the kaon mass 
(A+7r~7r + ). The last data sample would reveal excited 
A+ decays in which the tt~ has been misidentified as 
a AT~; no evidence for such decays is found. We also 
determined that misidentifed excited A+ decays result in 
peaks that are asymmetric in mass and are much broader 
than the signals observed in data. The wrong-sign data 
sample could show evidence of peaking backgrounds due 
to E° decays that are reconstructed with an additional 
kaon, but no such evidence is found. The A+ mass- 
sideband data sample and the generic-cc MC sample are 
searched for unexpected sources of peaking background. 
No sources of peaking background are found. We con- 
clude from these studies that the signals for S c (2980) + , 
£ c (3055)+, S c (3077)+, and S c (3123)+ arc not peaking 
backgrounds. These studies also confirm that the thresh- 
old function in Eq. (4) accurately models non-peaking 
background shapes. 



A+K^n- Results 



Projections of the A~^K®ir~ invariant-mass distribu- 
tion, for all five reconstructed A+ decay modes combined, 
are shown in Figs. 2(b) and (c). Projections of the com- 
bined compenents of the fit are also shown. The data and 
fit projections correspond to A'/(A+7r~) ranges within 
3.0 natural widths of the S c (2455)° mass (Fig. 2(b))and 
2.0 natural widths of the E c (2520)° mass (Fig. 2(c)) [2]. 
These M(A+7r~) ranges are delineated by the horizon- 
tal lines in Fig. 2(a). The fit shown in Fig. 2 accounts 
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M(A:K"7[ + ) (GeV/c 2 ) 



FIG. 1: (a) Scatter plot of M(A+tt+) versus M{A+K--rr + ), 
and (b,c) projections of the M(A+ K~n + ) distributions, in 
two ranges of M(A+7r + ), for data (points with error bars) 
summed over all five reconstructed Aj decay modes combined. 
The curves correspond to the two-dimensional fit results. The 
M(A^7r + ) ranges are (b) within 3.0 natural widths of the 
E c (2455) ++ mass, and (c) within 2.0 natural widths of the 
S c (2520) ++ mass. These M(A+7r + ) ranges are delineated by 
the horizontal lines on the scatter plot in (a). The projections 
of the combined background PDF components are illustrated 
by the dashed curves. 



for two excited charm-strange baryons: S c (2980)° and 
S c (3077)°. A two-dimensional binned x 2 probability of 
62% is calculated (in the same way as described for the fit 
to A+K-TT+ data) for the fit when both the S c (2980)° 
and S c (3077)° PDF components arc included, while a 
54% x 2 probability is calculated when only the S c (3077)° 
PDF component is included. Estimated statistical sig- 



nificances of 1.7 a (Aln£ = 8) and 4.5 a (Aln£ = 21) 
are found for the S c (2980)° and S c (3077)° states, re- 
spectively, based on changes in ln-likclihood values and 
the number of free parameters in each signal PDF com- 
ponent. Despite the low statistical significance of the 
S c (2980)° signal, we use the prior assumption that this 
state exists as the isospin partner of the S c (2980) + and 
include the S c (2980)° PDF components for all measure- 
ments based on the sample of A+K®ir~ candidates. 

The S c (2980)° PDF includes both a nonresonant 



A+K\ 



component and an intermediate-resonant 



£ c (2455)°if° component. The 2 C (3077)° PDF includes 
a nonresonant A+K®ir~ component, and intermediate- 
resonant S c (2455)°^° and ^(2520)°^° components; 
(78 ± 21 ± 5)% of the S c (3077)° signal is found 
to decay through the intermediate resonances. Of 
this intermediate-resonant signal, (44 ± 12 ± 7)% is 
£ c (2455) ++ -fr'~. No statistically significant signals are 
found for S c (3055)° or S c (3123)° states (the neutral 
isospin partners of S c (3055) + and S c (3123) + ). 



C. Systematic Uncertainties 

Several sources of systematic uncertainty are investi- 
gated and quantified for the measurements of charm- 
strange baryons decaying to A+A"~7r + and A+K®ir~ final 
states; they are listed in Table II. Uncertainties in PDF 
modeling are estimated from additional fits to the data 
with different mass resolution parameters, background 
shapes, and phase-space threshold masses. Because of 
their relatively low statistical significances, we include 
the effects of excluding a S c (2980)° or S c (3123)+ PDF 
component in the relevant fits. Systematic uncertainties 
on the measured mass, associated with SVT misalign- 
ments, angular dependence of tracking performance, en- 
ergy loss in detector material, the magnetic field, and ma- 
terial magnetization, were extensively studied in BABAR 
for a precision measurement of the A+ mass [15]. Since 
the Q- values of the excited charm-strange baryon decays 
in this analysis are similar to those of the A+ decays used 
in the mass measurement, we assign the same systematic 
error of 0.14MeV/c 2 . This uncertainty is much smaller 
than those from other sources. 

The systematic uncertainties are added in quadrature 
and are listed along with the measurements of masses, 
widths, yields, and intermediate-resonant decay fractions 
in Table I. To determine the effect of systematic uncer- 
tainties on the estimated significance of the S c (3055) + 
and S c (3123) + states, the change in ln-likelihood is recal- 
culated with the same PDF modeling changes described 
in the previous paragraph. Only for the 5 C (3123) + do 
any of these changes reduce the significance of the signal; 
the modified background shapes reduce the significance 
from 3.6 to 3.0 a (Aln£ = 12). This lower significance 
is listed in parentheses in Table I. 
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M(A C + K s 7i") (GeV/c 2 ) 



FIG. 2: (a) Scatter plot of M(A+tt _ ) versus M(A+K^), and 
(b,c) projections of the M(Aj K°-k~) distribution, in two ranges 
of M(A+7r~), for data (points with error bars) summed over 
all five reconstructed A J decay modes combined. The curves 
correspond to the two-dimensional fit results. The M(A^n~) 
ranges are (b) within 3.0 natural widths of the E c (2455)° mass, 
and (c) within 2.0 natural widths of the E c (2520)° mass. These 
Af (A+7T - ) ranges are delineated by the horizontal lines on the 
scatter plot in (a). The projections of the combined background 
PDF components are illustrated by the dashed curves. 



D. Products of Cross-Section and Branching 
Fractions 

Reconstruction efficiencies are estimated from JET- 
SET74 simulations of the e + e~ production of excited 
charm-strange baryons. The efficiencies arc studied 
as a function of the kinematic variables M(A^K) 2 , 
M(iC7r ± ) 2 , and A/(A ( f7r ± ) 2 , where K denotes K~ or K° s . 



TABLE I: Measured masses, widths, yields, resonant de- 
cay fractions, and significances for E c (2980) + , H c (3055) + , 
H c (3077) + , and H c (3123) + baryons decaying to A+K~tt + , 
and H c (2980)° and S c (3077)° baryons decaying to A+K^n~. 
The first errors are statistical and the second errors are sys- 
tematic. All signal candidates are required to have p* greater 
than 2.9GeV/c. The rows labeled "Resonant" give the frac- 
tion of each signal that decays through intermediate reso- 
nances. The E c (3077) + resonant fraction is a 90% confidence- 
level lower limit. The row labeled "E c (2455)" gives the frac- 
tion of resonant decays that proceed through E C (2455)AT . The 
only state for which systematic uncertainties reduce the cal- 
culated significance at the stated precision is the H c (3123) + ; 
the lower significance is given in parentheses. 





E c (2980)+ 


H c (2980) u 


Mass (MeV/c^) 


2969.3 ± 2.2 ± 1.7 


2972.9 ±4.4 ±1.6 


Width ( MeV) 


27±8±2 


31±7±8 


Yield 


756 ± 178 ± 104 


67 ± 33 ± 29 


Resonant (%) 


55±7±13 




Significance 


> 9.0 a 


1.7a 




E c (3077)+ 


S c (3077) u 


Mass (MeV/c*) 


3077.0 ± 0.4 ± 0.2 


3079.3 ±1.1 ±0.2 


Width ( MeV) 


5.5 ±1.3 ±0.6 


5.9 ±2.3 ±1.5 


Yield 


403 ± 54 ± 27 


90 ± 22 ± 15 


Resonant (%) 


> 80 


78 ± 21 ± 5 


E c (2455) (%) 


45 ± 5 ± 5 


44 ± 12 ± 7 


Significance 


> 9.0a 


4.5 a 




E c (3055)+ 


H c (3123)+ 


Mass (MeV/c^) 


3054.2 ± 1.2 ± 0.5 


3122.9 ±1.3 ±0.3 


Width ( MeV) 


17±6±11 


4.4 ± 3.4 ± 1.7 


Yield 


218 ± 53 ± 79 


101 ± 34 ± 9 


Significance 


6.4a 


3.6 a (3.0 a) 



A slight 6% relative variation in efficiency is found only 
at large values of Af(A+7r + ) 2 . Therefore, we use an ef- 
ficiency for each A+K~it + candidate that depends on 
the measured value of M(A+7r + ) 2 . Reconstruction ef- 
ficiencies for states that decay through A+K~ir + range 
from about 1% to 9%, depending on the reconstructed 
A+ decay mode. Reconstruction efficiencies for states 
that decay through A+/v!?7r _ range from about 0.4% to 
3.2%. 

For each excited charm-strange baryon state, the num- 
ber of baryons produced in the BABAR detector that de- 
cay to A+i\"~7r + or A+A'g7r~, where A+ — > pK~ir + , is 
estimated using the signal yields for all five A c decay 
modes, the reconstruction efficiencies, the A+ branching- 
fraction ratios [2], and a "Best Linear Unbiased Esti- 
mate" (BLUE) method [16]. The BLUE method ac- 
counts for correlated errors between the five signal yields 
from the five reconstructed A+ decay modes, as well as 
correlated errors between the five estimated efficiencies. 
The estimated number of excited charm-strange baryons 
produced with p* > 2.9GeV/c and the integrated lumi- 
nosity are used to calculate the product of cross-section 
and branching fractions for each excited charm-strange 
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TABLE II: Systematic uncertainties on the masses, widths, 
yields, resonant fractions (-Ri), and E c (2455) resonant frac- 
tions (R2) determined from the samples of A^K~n + and 
A^K°tt~ candidates. Uncertainties associated with mass res- 
olution, background-PDF shapes, phase-space thresholds, the 
inclusion of the E c (3123) + or H c (2980)° signal shapes, and the 
mass scale are listed. The systematic errors from each source 
are added in quadrature. NA indicates that a source of sys- 
tematic uncertainty is not applicable. 





Mass 


Width 


Yield 


Ri 


#2 




(MeV/c 2 ) 


(MeV) 


(%) 


(%) 


(%) 














Mass Resolution 


1 P. f: 

ztU.u 


-4-1 Q 

ztl.o 


1 A 

zt 4 


1 A 

zt 4 


IN A 


Background Shape 


-4-i 

zt l.U 


ztU.o 


4- Q 
zt 


4- S 
zt 


IN A 


Phase-Space Thresh. 


zt l.z 


in 
ztU.o 


4-1 O 
ztlU 


1 0, 

zt y 


M A 
IN A 


S!i PTifil Tnpln Qinn 


±0.4 


±0.4 


± 3 


± 3 


NA 


]VIass Scale 


±0.1 


NA 


NA 


NA 


NA 


Total 


±1.7 


±1.5 


±14 


±13 


NA 


ZL c {dUOQ) 












Mass Resolution 


in a 
ztU.4 


4-7 Q 


zt ZD 


NT A 
IN A 


M A 
IN A 


Background Shape 


in 
ztU.o 


1 a 
ztO.o 


4-00 

ztZz 


NT A 


M A 
IN A 


PVifl^p-Slnfipp TVirptjVi 

± ±±CLo^ kJ V-f CLVj^Z _1_ ILL CDll . 


±0.1 


±4.2 


±11 


NA 


NA 


Signal Inclusion 


±0.0 


±1.4 


± 5 


NA 


NA 


Mass Scale 


±0.1 


NA 


NA 


NA 


NA 




-4-0 ^ 

Z1ZU.O 


4-1 1 


ZIZOU 


N A 


MA 
IN /\ 


[ou I i ) 












Mass Resolution 


in 1 -1 
ztU. 1 1 


in a 
ztU.4 


1 O O 

zto.Z 


zto.U 


4-0 

ztz.O 


R fi ok crvdi lTiri Shanp 


±0.10 


±0.2 


±0.1 


±1.2 


±2.7 


Phase-Space Approx. 


±0.09 


±0.3 


±3.0 


±1.4 


±1.9 


Signal Inclusion 


±0.05 


±0.2 


±5.0 


±4.6 


±2.4 




±0.14 


NA 


NA 


NA 


NA 


Total 




4-0 

ZIZU.U 




4-fi 1 


±4.8 


l^c l^olZOj 












Mass Resolution 


in 
ztU.o 


4-1 K 
±1.0 


1 cr n 
zto.U 


TVT A 
IN A 


M A 
IN A 


R fi pW ornii n rl SlVifiT^P 


±0.2 


±0.6 


±6.9 


NA 


NA 


Phase- Space Thresh. 


in 1 
ztU.l 


in tr 
ztU.O 


1 n 
zto.U 


v \ 

IN A 


M A 
IN A 


Mass Scale 


±0.1 


NA 


NA 


NA 


NA 


Total 


±0.3 


±1.7 


±8.9 


NA 


NA 


H c (2980) u 












Mass Resolution 


±0.6 


±7.1 


±16 


± 1 


NA 


Background Shape 


±1.3 


±3.7 


±37 


±14 


NA 


Phase-Space Thresh. 


±0.7 


±1.7 


±12 


±17 


NA 


Mass Scale 


±0.1 


NA 


NA 


NA 


NA 


Total 


±1.6 


±8.2 


±42 


±22 


NA 


H c (3077) u 












Mass Resolution 


±0.01 


±0.3 


± 1 


±0.4 


±0.3 


Background Shape 


±0.12 


±0.2 


± 1 


±4.1 


±3.1 


Phase-Space Thresh. 


±0.03 


±0.1 


± 1 


±2.4 


±0.2 


Signal Inclusion 


±0.02 


±1.5 


±17 


±0.9 


±6.3 


Mass Scale 


±0.14 


NA 


NA 


NA 


NA 


Total 


±0.19 


±1.5 


±17 


±4.9 


±7.0 



baryon: 

a(e+e- - E* C X)B(E* C - F)B(A+ -> p^-7T+) , (5) 

where S* is an excited charm-strange baryon, A" is the 
unreconstructed portion of the e + e~ event, and Y is 
A+K~w + or A+AT°7r _ . With these calculations, it is as- 
sumed that 34.6% of the K° mesons decay as K® — ► 



TABLE III: The products of cross-section and branching frac- 
tions (Eq. (5)) for excited charm-strange baryons produced 
with p* > 2.9 GeV/c and decaying to three-body final states. 
The first errors are statistical and the second errors are sys- 
tematic. Upper limits are at 90% confidence level. 



1 — 'c 




A+A u tt- 


H c (2980) 


(11.8 ± 3.4 ± 2.2) fb 


< 15 fb 


H c (3055) 


(2.2±1.2±0.7)fb 


< 7fb 


H c (3077) 


(8.1 ± 1.2±0.8)fb 


(6.2±2.1±1.5)fb 


H c (3123) 


(1.6±0.6±0.2)fb 


< 1.4 fb 



7r + 7r [2]. The measured products of cross-section and 
branching fractions are listed in Table III. 

Upper limits at 90% confidence level are calculated on 
the product of cross-section and branching fractions for 
the neutral states S c (2980)°, S c (3055)°, and H c (3123)°. 
The upper limits are determined from numerical inte- 
gration of posterior probability densities calculated from 
the likelihood distributions, with a uniform positive prior 
for the product of cross-section and branching fractions 
greater than zero and a zero prior below. Gaussian con- 
straints on efficiencies and A+ branching-fraction ratios 
are included in the likelihood functions; these modified 
likelihood functions constrain the ratios of fitted yields 
from each A+ decay mode and incorporate systematic 
uncertainties from these ratios of yields. Also, in calcu- 
lating each upper limit, the signal shape parameters of 
the baryon are given Gaussian constraints based on the 
measured values in Table I. For S c (3055)° and S c (3123)° 
parameters, we use the values of the parameters mea- 
sured for their isospin partners. All Gaussian constraints 
involve a factor analogous to that in Eq. (3). The up- 
per limits on the products of cross-section and branching 
fractions are listed in Table III. 



IV. DECAYS TO A+K% AND A+K~ 

The studies of decays to the two-body final states 
A^K° S and A+A'~ are based on one-dimensional fits to 
the distributions of M(A+if°) and M(A+K~). In each 
case, the invariant-mass range between 2.91GcV/c 2 to 
3.15GcV/c 2 is fit. The PDF component that accounts 
for a possible resonance (P s ) is a double- Voigtian reso- 
nance shape with free parameters a, which are the mass 
and natural width of the signal. In the likelihood func- 
tion, the means and widths of the resonance shapes are 
given Gaussian constraints (Eq. (3)) based on the values 
measured with three-body decays. For each of the five 
A+ decay modes, the background distribution in A+if° is 
modeled with a first-order polynomial (P&). Background 
distributions in A+K~ are modeled with functions pro- 
portional to Eq. (4). 

Invariant-mass resolutions and reconstruction efficien- 
cies are calculated from MC samples. Mass resolutions 
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FIG. 3: Distributions of (a) M(A+Jfg) and (b) M(A+K~), for 
all five reconstructed A+ decay modes combined. The curves 
correspond to fits with no signal PDFs. Signals were searched 
for around the invariant-mass values indicated by the dashed 
vertical lines. 



are about 2.5MeV/c 2 . Reconstruction efficiencies range 
from about 0.7% to 5.0% for signals that decay through 
A+ifg, depending on the reconstructed A+ decay mode, 
and from about 1.5% to 12% for signals that decay 
through A+A'~. 

Neither the A + AT° nor the A+A'~ invariant-mass dis- 
tributions exhibit evidence for statistically significant 
peaking structures corresponding to any of the charm- 
strange baryons found in the three-body final states. Fig- 
ure 3 shows background-only fits to the A+K° and A + AT - 
invariant-mass distributions for all five reconstructed A+ 
decay modes combined. The same method described 
for three-body final states is used to determine 90% 
confidence-level upper limits for the products of cross- 
section and branching fractions described by Eq. (5), 
where Y is now A+A' or A + A'~. The upper limits on 
the products of cross-section and branching fractions are 
listed in Table IV. 



V. DECAYS TO A+K°n-ir+ AND A+K~iv + it~ 

The studies of decays to the four-body final states 
A+ATg7r~7r + and A+A'~7r + 7r~ are based on fits to 
the one-dimensional distributions of M(A+A"g7r~7r + ) 



TABLE IV: 90% confidence-level upper limits on the products 
of cross-section and branching fractions (Eq. (5)) for excited 
charm-strange baryons produced with p* > 2.9 GeV/c and de- 
caying to two-body and four-body final states. NA indicates 
decays that are not kinematically allowed. NO indicates that 
corresponding three-body decays are not observed. The four- 
body final-state upper limits only use integrated yields up to 
about 30 MeV/c 2 above kinematic threshold. 



1 — 'c ^ 


A+A* u 


A+A- 


A+A" U 7r-7r+ 


A+A' 7T + 7T 


E c (2980) 


< 12 fb 


< 10 fb 


NA 


NA 


H c (3055) 


< 9fb 


NO 


NA 


NA 


H c (3077) 


< 2.9fb 


< 1.2 fb 


< 0.4fb 


< 0.1 fb 


E c (3123) 


< 2.7fb 


NO 


< 1.4 fb 


NO 



and M(A + A' 7r + 7r ). The invariant-mass region from 
the minimum kinematic threshold (« 3063MeV/c 2 for 
A+A:07r-7r+ and « 3059MeV/c 2 for A+K-jt+tt-) to 
250MeV/c 2 above the threshold is fit in each case. The 
PDF component that accounts for possible resonances 
is a double- Voigtian resonant function multiplied by a 
third-order polynomial that models the rapidly chang- 
ing four-body phase space. In the likelihood function, 
the means and widths of the resonance shapes are given 
Gaussian constraints (Eq. (3)) based on the values mea- 
sured with three-body decays. The rapidly rising four- 
body phase-space function and the proximity of any res- 
onances to the kinematic thresholds lead to possibly sig- 
nificant contributions of S c (3077)° and S c (3077) + signals 
throughout each fitted 250 MeV/c 2 range. In quoting up- 
per limits for S c (3077)° and S c (3077) + , we consider the 
integrated yield up to 3093MeV/c 2 and 3089MeV/c 2 , re- 
spectively (« 30MeV/c 2 above threshold in each case). 
Each range includes the peaking signal but not the higher 
mass range where efficiencies and possible effects from 
intermediate-resonant decays are unknown. Background 
distributions are modeled by (M — a)(M — t) 2 , where M 
is the mass variable, t is the kinematic threshold, and a 
is a free parameter. 

Invariant-mass resolutions and reconstruction efficien- 
cies are again calculated from MC samples. Mass res- 
olutions are about 1.0MeV/c 2 . Reconstruction efficien- 
cies range from about 0.3% to 2.4% for signals that 
decay through A+A"g7r~7r + , depending on the recon- 
structed A+ decay mode. Reconstruction efficiencies 
range from about 0.8% to 5.5% for signals that decay 
through A^K~tt + tt~ . 

Neither the A+K®ir~ir + nor the A + AT~7r~7r + 
invariant-mass distributions exhibit evidence for statisti- 
cally significant peaking structures corresponding to any 
of the charm-strange baryons found in the three-body 
final states. Figure 4 shows background-only fits to these 
invariant-mass distributions for all five reconstructed 
A+ decay modes combined. The same method described 
for three-body final states is used to determine 90% 
confidence-level upper limits for the products of cross- 
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FIG. 4: Distributions of (a) M(A+^7r"7r + ) and (b) 
M(A^ K~ -k + tt~), for all five reconstructed decay modes 
combined. The curves correspond to fits with no signal PDFs. 
Signals were searched for around the invariant-mass values in- 
dicated by the dashed vertical lines. 



section and branching fractions described by Eq. (5), 
where Y is now A+A'°7r~7r + or A+A'~7r + 7r~. The upper 
limits on the products of cross-section and branching 
fractions are listed in Table IV. 



VI. SUMMARY AND CONCLUSIONS 

Invariant-mass distributions for six different final 
states are studied for evidence of excited charm-strange 
baryons. Four statistically significant signals are found 
for excited charm-strange baryons decaying to the final 
state A+AT-7T+. They are named S c (2980)+ (> 9er), 
S c (3055)+ (6.4(7), S c (3077)+ (>9a), and S c (3123)+ 
(3.8 cr), and their masses, widths, and the products 
of cross-section and branching fractions are measured. 
Intermediate-resonant E c (2455) ++ A'~ and nonrcsonant 
A+A"~7r + decays of the S c (2980) + are observed, while 
the !H C (3077) + is observed to decay mostly through 
the £ c (2455) ++ A:- and £ C (2520)++A"- intermediate- 
resonant states. The S c (3055) + and S c (3123) + signals 
are observed only in S c (2455) ++ A"" and E c (2520) ++ A'~ 
intermediate-resonant decays, respectively. For the 
A+A"s7r~ final state, the only statistically significant 
signal corresponds to the S c (3077)° baryon, which is 



presumably the isospin partner to the S c (3077) + ; its 
mass, width, and product of cross-section and branch- 
ing fractions are measured. Like the S c (3077) + state, 
the 5 C (3077)° is observed to decay mostly through the 
E C (2455)°AT° and S C (2520) AT° intermediate-resonant 
states. A 1.7 a enhancement in the final state A+AT"7r~ 
is measured to have a mass and width consistent with 
its being the neutral isospin partner of the S c (2980) + 
baryon. For the three-body final state A+AT!?7r - , 90% 
confidence-level upper limits are determined for the 
products of cross-section and branching fractions of 
S c (2980)°, S c (3055)°, and S c (3123)°. No statistically 
significant signals are found in the two-body and four- 
body final states. For states that are observed in three- 
body decays, we report 90% confidence-level upper lim- 
its for the product of cross-section and branching frac- 
tions for kinematically allowed two-body and four-body 
decays. 

The measured resonant and nonresonant decay rates 
and decay modes of these new excited charm-strange 
baryons might provide information on the internal quark 
dynamics. In both the resonant and nonresonant strong 
decays, the strange quark is contained in a kaon, sep- 
arate from the charmed baryon. Previously known ex- 
cited S c baryons have been observed only in decays to 
a lower-mass S c baryon plus a pion or photon. Theo- 
retically, there are several excited charm-strange baryon 
states with various spin, angular momentum, and radial 
excitation configurations [6-10]. These different theoret- 
ical states may offer explanations for differences in decay 
rates and decay modes, but the current theoretical and 
experimental information is not definitive enough to as- 
sign quantum numbers to any of the new excited charm- 
strange states. Future experimental studies of these and 
other states and decay modes, and further theoretical 
work, will help to clarify the properties of these new 
charm-strange baryons. 
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